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22 \ ABSTRACT 

Kf" ' We investigate the neighborhood of Topological Lattice Field Theories (TLFTs) 

■^ \ 
0\ , in the parameter space of general lattice field theories in dimension D >2, and dis- 

r^ I cuss the phase structures associated to them. We first define a volume-dependent 

"^ . TLFT, and discuss its decomposition to a direct sum of irreducible TLFTs, which 

_^ ' cannot be decomposed anymore. Using this decomposed form, we discuss phase 

structures and renormalization group flows of volume-dependent TLFTs. We find 
that TLFTs are on multiple first order phase transition points as well as on fixed 
>• . points of the flow. The phase structures are controlled by the physical states on 

/\ ' (-D — l)-sphere of TLFTs. The flow agrees with the Nienhuis-Nauenberg criterion. 

We also discuss the neighborhood of a TLFT in general directions by a perturba- 
tive method, so-called cluster expansion. We investigate especially the Zp analogue 
of the Turaev-Viro model, and find that the TLFT is in general on a higher order 
discrete phase transition point. The phase structures depend on the topology of 
the base manifold and are controlled by the physical states on topologically non- 
trivial surfaces. We also discuss the correlation lengths of local fluctuations, and 
find long-range modes propagating along topological defects. Thus various discrete 
phase transitions are associated to TLFTs. 
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1. Introduction 

The topological field theory (TFT) ' has been an important tool in theoreti- 
cal physics and pure mathematics. In theoretical physics it has been appearing in 
the investigations of the non-perturbative aspects of string theories and quantum 
gravities. The topological phase of two-dimensional quantum gravity was discussed 
by Witten . In three dimensions, he pointed out that Chern-Simons theories with 
certain gauge groups are interesting models for quantum gravity . The three- 
dimensional Chern-Simons theories have been also useful in two-dimensional con- 
formal field theories ^ and integrable models . TFTs are also interesting from 
mathematical point of view concerning topological invariants of manifolds ' . 

Recently lattice constructions of TFTs have been being discussed. They are 
simpler than the usual continuum formulations of TFTs in the sense that they are 
free of infinitely many local spurious degrees of freedom and divergences. This 
simpleness is expected to facilitate the systematic constructions and classifications 
of TFTs. 

In fact, in two dimensions, topological lattice field theories (TLFTs) have been 
constructed and classified extensively . It was shown that the class of TLFTs 
with degrees of freedom on links has one-to-one correspondence to semi-simple 
associative algebras ^ . The relation between this class of TLFTs and known 
topological matter field theories in two dimensions ' was discussed ' . In general, 
it was shown that, starting with Atiyah's axioms , any unitary TFT can be 
decomposed as a direct sum of irreducible TFTs, which can be realized as TLFTs 
with degrees of freedom on vertices. Here the notion of irreducible TFTs was 
first introduced, and they are TFTs which can not be decomposed anymore. This 
notion will play an important role in the present paper. 

In three dimensions, some systematic constructions of some series of TLFTs 
have been done, but any systematic classifications have not been done yet since 
they must be tightly related to the classification of three-dimensional manifolds. 
The first TLFT was constructed by Ponzano and Regge , motivated by the three- 



dimensional Regge calculus . The degrees of freedom of their TLFT are the la- 
bels of the irreducible representations of SU{2). Their TLFT was generalized to 
a series of models with the quantum groups SUq{2) by Turaev and Viro . Di- 

[221 

jkgraaf and Witten constructed models with discrete gauge groups. A relation 
between the Turaev- Viro model and 2D conformal field theories was pointed out, 
and some generalizations of the model were done ' . Boulatov ' constructed 
tensor-model versions of the Turaev- Viro model. (See also [26]-[31] for other dis- 
cussions concerning Turaev- Viro model.) The relations to continuum formulations 
have not been studied enough. But it was shown by Turaev that the Turaev- Viro 
model is equivalent to the square of the Chern-Simons theory with the gauge group 
SU{2). Ooguri and the present author ' discussed the equivalence between the 
Ponzano- Regge model and the Chern-Simons theory with gauge group /SO (3). A 
new class of TLFTs was constructed by Chung, Fukuma and Shapere . This 
class of TLFTs is defined on arbitrary 3D cell complexes, not only on tetrahedral 
lattices. These TLFTs were shown to have one-to-one correspondence with a class 
of Hopf algebras. 

In four dimensions, Ooguri "^ generalized the Boulatov's argument to four di- 
mensions and obtained the 4D analogies of the Ponzano- Regge- Turaev- Viro model. 
He argued that the 4D analogue of the Ponzano- Regge model is the lattice realiza- 
tion of the i?F-theory 

The recent developments on TLFTs have been concentrated only on construc- 
tions and classifications of TLFTs, and seem to uncover physical content. Since 
TLFTs are special cases of general lattice field theories, the investigation of the 
positions of TLFTs in the parameter space of general lattice field theories will 
clarify the specialties of TFTs and the background reason of their appearances 
in the investigations of quantum gravity and string theory. The present author 
investigated the phase structure around 2D TLFTs by a perturbative scheme and 
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found that they are in general on multiple first order phase transition points 
A usual continuum field theory describes the long-range modes near a fixed point 
of a renormalization group flow on a continuous phase transition surface. Thus, 



complimentary to usual continuum field theories, TLFTs(TFTs) might be useful 
in describing the physics near a fixed point on a discrete phase transition surface. 

This paper is an extension of the present author's previous paper to dimension 
D > 2. The paper is organized as follows. 

In section 2, we will define volume-dependent TLFTs in any dimension D{> 2), 
and discuss their decompositions under certain conditions. We will review the 
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physical structures of TLFTs . Then we will decompose a volume-dependent 
TLFT satisfying certain conditions into a direct sum of irreducible TLFTs with 
volume-dependent numerical factors. Using this decomposed form, we will show 
that a TLFT is in general on a multiple first order phase transition point, which 
is also a fixed point of the renormalization group fiow in the parameter space of 
volume-dependent TLFTs. As an example, the procedure of the decomposition is 
shown in the two-dimensional lattice QCD ' ' . 

In section 3, we will investigate the neighborhood of a TLFT by a perturbative 
method. We will discuss the method in detail, which is a kind of cluster expansion. 
Using this method, we will investigate the phase structure near a TLFT in general 
directions, and again will find the structure of first order phase transitions in the 
first order of the cluster expansion. The number of the different phases around 
it is shown to be equal to the number of the irreducible TLFTs of which the 
TLFT is a direct sum. The higher order contributions do not change this structure 
qualitatively, provided the perturbative treatment is valid around a TLFT and 
that the base manifold is topologically trivial. As a simple example, we will discuss 
the neighborhoods of the TLFTs in the parameter space of the 2D Q-state Potts 

1 ,[13] 

model 

In section 4, we will investigate the phase structure near the discrete group Zp 
analogue ' '' of the Turaev-Viro model on a topologically non-trivial manifold in 
three dimensions. This model is an irreducible TLFT, hence is not on a first order 
phase transition point. But we will find that the model is in general on a higher 
order discrete phase transition point, which is controlled by the physical states 



on topologically non-trivial surface. To show it, the higher orders of the cluster 
expansion are essential. 

In section 5, we will give summary, comments and discussions. 



2. Volume-dependent topological lattice field theories 

For later discussions, it is more convenient to enlarge the parameter space of 
TLFTs. We will begin with defining volume-dependent TLFTs, which do not have 
any local degrees of freedom like TLFTs but are dependent on the volume of the 
manifold. 

The rough description of the usual constructions of TLFTs is as follows. Topo- 
logical field theories are field theories which are invariant under local deformations 
of background metric. Physical quantities of such theories depend only on the 
topology of the underlying manifold. In TLFTs, the independence of the back- 
ground metric of TFTs corresponds to the independence of the lattice structure. 
Consider a simplicial decomposition of a given manifold. Usually lattice theories 
are defined by associating degrees of freedom to the simplices and assigning local 
weights. The partition function is defined as the summation of the products of 
the local weights over the degrees of freedom. The partition function of a TLFT 
must be invariant under the change of the simplicial decomposition of the manifold. 
The simplicial decompositions of a manifold can be generated by a certain set of 
local moves in any dimension ' . These local moves constrain the local weights. 
Solving these constraints, one obtains TLFTs. 

In the present paper, we are not interested in concrete models of TLFTs, but in 
the general structures of TLFTs in an arbitrary dimension. Thus, for convenience, 
we use a more formal definition of TLFTs. Following Atiyah , we define a TLFT 
as a set of partition functions associated to simplicially decomposed manifolds 
possibly with boundaries. The partition functions should be independent of the 
simplicial decompositions of the insides of the manifolds. Hence these partition 



functions must satisfy certain consistency conditions under gluing operations of 
the manifolds. 

2.1 Definition of volume- dependent TLFTs 

The lattices we consider in this paper are those obtained by simplicial decom- 
positions of D-dimensional oriented connected manifolds possibly with boundaries. 
The lattices of boundaries can be colored by associating to each /-dimensional sim- 
plex (/ = 0, 1, ■ ■ ■ , D — 1) an element xi of a finite index set Xi. 

The ingredients of a volume-dependent TLFT are complex-valued partition 
functions. A partition function is associated to a manifold M which is decomposed 
with A D-dimensional simplices and possibly has colored boundaries S(c); 

Z(M;A;S(c)). (2.1) 

Here c denotes symbolically the coloring. The partition function is independent 
of the way how the inside of the manifold is decomposed with A D-dimensional 
simplices but depends on the number A and the topology of M as well as the 
simplicial decomposition and the coloring of the boundaries. 

The partition functions (2.1) must satisfy constraints under gluing operations 
of lattices. Let (M; A; S) denote a lattice with uncolored boundaries. Consider a 
lattice (M3; A3; S3) obtained by identifying pairwise some of the {D — l)-simplices 
of the boundaries of (Mi; Ai; Si) and (M2; A2; S2) (Fig.l). Let b denote the sub- 
boundaries identified pairwise. Then A^ = Ai + A2, M3 = Mi U^ M2 and S3 = 
Si Uf, S2, symbolically. We demand that the partition functions for the three 
lattices should satisfy the following constraint under the gluing operation: 

Z{Ms- As; S3(C3)) = J2 ^(^15 Ar, Si(ci))(7^^"^^^Z(M2; A2; S2(c2)). (2.2) 

Here 5'^"^^'' is the gluing tensor, and Cii,{i = 1,2) are the colorings associated to the 



identified simplices of tlie sub-boundaries . Tlie summation is over tlie colorings 
associated to tlie simplices which become inner simplices after the gluing operation. 

We demand also that the partition functions must satisfy the following con- 
straint under a self-gluing operation: 

Z(M2; A- S2(C2)) = Y, C^'ZiMi; A; Si(ci)). (2.3) 

C-inn 

Here (M2; A; S2) is obtained by identifying some of the {D — l)-simplices of the 
boundaries Si (Fig. 2). 

2.2 Laplace decomposition 

In this subsection, we will decompose a volume-dependent TLFT satisfying 
a certain condition into a direct sum of TLFTs by applying an analogue of the 
Laplace transformation to the volume-dependence of the partition functions. Then, 
in the following subsections, we will decompose each TLFT into a direct sum of 
irreducible TLFTs. 

Let T denote a D-dimensional simplex, and dT its boundary lattice. We 
consider only a TLFT satisfying the following condition: A discrete analogue of 
the Laplace transformation is applicable to the volume dependence of the partition 
function for a D- dimensional ball B^ with dT^ as its boundary. 

ZiB^"; A; dT^'ic)) = ^ Z,(i?^; 9r^(c))a;/, (2.4) 

i 

where Ui 's are complex numbers and i is the label to distinguish different Ui 's. 

First we will prove that the Zi(B ,dT (c))^s are orthogonal. Consider two 
lattices {B , Ai; dT ) {i = 1, 2) and a lattice (B ; A\ ^ A% S) obtained by gluing 
the two lattices identifying one of the [D — l)-simplices on dT of each lattice 



• In most of concrete models of TLFTs, the gluing tensor is the product of the Kronecker's 
deltas of the colors of the identified simplices with some numerical factors. 



(Fig. 3). Since the partition functions must satisfy the constraint (2.2) for the 
gluing operation, one obtains 

Z(i?^; Ai + A2; S(C3)) = J2 ^(^""5 Ai; dT^{ci))g'^^,'^lZ{B^- A2; dT^{c2)). 

Cinn 

(2.5) 
Using the assumption (2.4), the right-hand side is rewritten as 

Y, E ^*(^''' dT''{ci))g''^'h'^\Z,{B''- 9T^(c2))c./^c./^ (2.6) 

Since the left-hand side of (2.5) is a function of Ai + A2, so must be (2.6). Thus 
the Zj's must be orthogonal, and we obtain the Laplace transformation for the 
left-hand side of (2.5) with 

(5y Z,(5^; S(C3)) = Y. Z,{B^; dr^{ci))g'^^,'^^Z,{B^; dr^{c2)). (2.7) 

Cinn 

In the self-gluing operation (2.3), if the Laplace transformation is applicable 
to the right-hand side, then so is the left-hand side. The Zj's in the both sides are 
related by 

Zi{M2;^2{c2)) = 5^(?^'"'"Z,(Mi;Ei(ci)). (2.8) 

Any partition function can be obtained by gluing operations of Z^-simplices and 
self-gluing operations. Thus repeating the same discussions as above, we obtain 
that the Laplace transformation using the same Ui^s is applicable to the volume- 
dependence of any partition function. Then the constraint (2.2) is rewritten using 
Zj-i s as 

<5yZ,(M3;S3(c3)) = 5^Z,(Mi;Si(ci))(?;;"^^''Z,(M2;S2(c2)). (2.9) 

C-inn 

The volume-dependence has disappeared in (2.8) and (2.9). Since Zj's are 
orthogonal for different i's, the constraints (2.8) and (2.9) are the defining con- 
straints for the partition functions of a TLFT for each i, as will be explained in 



the following subsection. Thus we have shown the following. A volume- dependent 

TLFT with (2.4) is a direct sum of TLFTs with some volume-dependent numerical 

factors: 

N 
Z{M; A; S(c)) = J^ Z,(M; S(c))cui^, (2.10) 

where N is a finite number , and Zi 's are the partition functions of TLFTs. 

2.3 Definitions concerning TLFTs 

To prepare to decompose (2.10) further, we will discuss some elementary defi- 
nitions concerning TLFTs in this subsection. 

The ingredients of a TLFT are complex- valued partition functions. A partition 
function of a TLFT is associated to a manifold M which is simplicially decomposed 
with some Z^-simplices and possibly has colored boundaries: 

Z(M;S(c)). (2.11) 

This partition function is independent of the way how the inside of M is decom- 
posed and must satisfy the constraints (2.8) and (2.9). 

The Hilbert space TYs associated to a simplicially decomposed boundary S is 
the module freely generated by all possible colorings of S over complex numbers. 
Thus a wave function in Tis is a complex- valued function of colorings on S;7is = 
{$e(ce)}- The Hilbert space is finite dimensional since the index set Xi is finite. 

Consider a, {D — 1) -dimensional closed manifold S, and M = S x [0, 1]. Take 
simplicial decompositions Si and Sg for the two boundaries of M so that one can 
take a simplicial decomposition of M consistent with Si and Sg. Here the symbol 



• The finiteness can be shown as follows. In the following subsection, we will discuss the 
physical Hilbert spaces of TLFTs and volume-dependent TLFTs. It will be shown that the 
dimension of the physical Hilbert space on the {D — l)-sphere is non-zero and finite. Thus 
N must be finite. 



* is for the reverse of orientation. We define a linear map which maps a state in 
Tisi ^o ^ state in 7ij]2- 



<^E2(c2)=i'E2Ei(c2,Ci)$Si(ci), 

(2.12) 
Ps,s,(c2,ci) =^^f Z(S X [0,l];S^(4),Ei(ci)), 



where the repeated indices for the colorings are summed over. From now on, it is 
supposed that the repeated indices for colorings are summed over unless otherwise 
stated. 

Since the manifold M = S x [0, 1] can be obtained by gluing two M's identifying 
the boundary S of one M and S* of another M, the linear maps satisfy the following 
equation: 

^S3Si(c3,Cl) = Ps3S2(c3,C2)Ps2Si(c2,Cl). (2.13) 

This equation (2.13) is derived from the definition (2.12) and the constraint (2.9) 
under the gluing operation above. Taking S3 = S2 = Si, we see that -PsiEi = 
-PeiEi-PsiSh ^-e., the map -PsiSi is a projection map. 

The physical Hilbert space 7Y|, ^^ is defined as the invariant subspace of Ti^i 
by the projection operator -PeiSi- The physical Hilbert space is finite dimensional 
since the Hilbert space is finite dimensional. Take a basis of the physical Hilbert 
space and its dual basis: 





0L 


.(^) 


e nl^y' 
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0Si 


^{C) 


e v$y' 


* 
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0Si 


.(c)0i 


M) 


= ^l 
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=1,---: 


AiH'ntr 



(2.14) 



Taking S3 = Si in (2.13), one finds -PeiSi = -PsiS2-Ps2Si- Thus -PS2S1 is a one-to- 
one map from 7Y|, ^^ to 7i|, ^^ . Using the one-to-one maps -PsiS2 and -Ps2Sn one 
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can obtain a natural basis and a natural dual basis for Tiv ^^'■ 

0E2(c2)=i'E2Si(c2,Ci)0i^^(ci), 

(2.15) 
0S2 i(c2) = 0S, j(ci)PsiS2(ci,C2). 

From now on, the physical states labeled with the same index are supposed to be 
related in this way. 

Consider a manifold M with boundaries S|, ■ ■ ■ , S* , S^+i, ■ ■ ■ , S^ with certain 
simplicial decompositions. Following (2.12), define 

PM;S*,-,E* ;S„+i,-,E„(ci, ■■■ ,Cn] Cn+1, " " " , Cm) 

2.16 

= C • ■ -C"^!^; ^^(4), ■ ■ ■ , s;(c;), s„+i(c„+i), • . ■ , s,„(cm)), 

which can be regarded as a linear map from ?is„+i <S> • • ■ <S> '^^Sjn to Tisi ® ■ ■ • ® 7is„ • 
Since gluing manifolds Sj x [0, 1] (? = 1, ■ ■ ■ , m) to M by identifying Sj's and S^'s 
does not change the manifold M, one can show, using the definitions (2.12) and 
(2.16) and the constraint (2.9) under the gluing operation, that 

-PM;E*,-,E*;S„+i,-,S™(ci, ■■■ .Cn] C„+i, ■ ■ ■ , Cm) 

= ^SiEi(ci, c'l) ■ ■ ■ Ps„S;(Cn, C^)PM;Er,-,E^*;E;+i,-,S;„(c'l, " " " , C^; cj^+i, ■ ■ ■ , cj^) 

X ^S;^+iS„+i(Cn+i, Cn+l) ■ " ■PE^E„(Cm, Cm), 

(2.17) 

where Sj and S^ denote the same boundary possibly with different simplicial de- 
compositions. The equation (2.17) is a generalization of (2.13). Thus, taking 
Si = S^ in (2.17), the linear map (2.16) is in fact a map among the physical 
Hilbert spaces; T^g^^'J ® ■ ■ ■ ® Tf^y' to T^g^^ ® ■ ■ ■ ® 7^g^^ This fact leads us 
naturally to the following definition of a physical correlation function: 

Z{M- (Si),„ ■ ■ ■ , (S„),„, (S,+i)^"+\ ■ ■ ■ , (Sm)*-) 

= 0Ei ii(ci) ■ ■ ■ 0S„ i„(Cn)PM;S*,-,S*;E„+i,-,E„(ci, " " " , Cn; Cn+1, " " " , Cm) 

(2.18) 
The remarkable property of a physical correlation function is that it is independent 
of the simplicial decompositions of its boundaries. This fact can be proven by 
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using the definition of tlie basis of the physical Hilbert space (2.15) and (2.17). In 
particular, for any S, 

ZiEx [0,1]; {E)„im = si, (2.19) 

because the projection map (2.12) is the identity on 7i|, ^'*. 

2.4 Decomposition of a TLFT 

Here we will discuss the decomposition of each TLFT in (2.10) into irreducible 
TLFTs. An irreducible TLFT is defined as a TLFT which cannot be decomposed 
further as a direct sum of TLFTs. This is an extension of the work by Durhuus and 
Jonsson , where they discussed the decomposition of a unitary two-dimensional 
TFT into a direct sum of irreducible TFTs. 

The discussions will be roughly as follows. First we will discuss the correlation 
functions of S s, and will find a convenient basis of the physical Hilbert space 
on S , which we call the canonical basis. This basis was first introduced by 
Verlinde in the discussions of fusion rules of 2D conformal field theories, and 
was applied to 2D TLFTs in [13] - [15] and [16]. Then we will show that any 
partition(correlation) functions of a TLFT can be expressed as a sum of the par- 
tition(correlation) functions labeled by the labels of the canonical basis and that 
these labeled partition functions are orthogonal. 

Here we begin with the notations to be used. Bq denotes S with a certain 
simplicial decomposition. Since, in the physical correlation functions, the simplicial 
decompositions of boundaries are irrelevant, we use 0* for the physical state 0^ 
and (pi for (pSo i- We also use (S)*(c) for (pj^{c) and (S)j(c) for 0|. ^(c), for short. 

Define the following symmetric tensors which reverse the orientations of the 
physical states on S^~^: 

Tjij = Z{TsD-i;(t)i,(t)j) 

= {BoUc,)g'j^^^Z{Tsn-v,B*o{c2),B*o{csW-^'''{Bo)j{c4), (2.20) 



V'^ = Z{Tsn-^; 



^y = 7(T^r. ,;0*,0J) 
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where Tgo-i = S ^ x [0, 1]. Since the projection operator (2.12) is the identity 
on the physical Hilbert space, the following equation holds for any S: 

5;](Sy(ci)(S),(c2)=PEs(ci,C2), (2.21) 

i 

From now on it is supposed that contracted indices labeling the physical states 
are summed over unless otherwise stated. Using (2.21), (2.19) and (2.9) for the 
operation that two T^-d-is are glued to make one Tgo-i, we obtain 



'^■^ = 5t (2.22) 



VijV' 



which shows that the tensors (2.20) are not degenerate. These tensors can be used 
to change the cases of the indices of (pi and 0*. The proof is as follows: 

Z(M;f ,S(c)) = Z(T5z,-i;0\5o(ci))(7g;^Z(M;i?o*(c2),S(c)) 

= ZiTsn-^;<p\ 5o(ci))Pb„Bo(ci, C2)g'i;;'Z{M; B*o{c,), S(c)) 

= ZiTso-v,<p\ Bo{ci))<P^ici)<Pjic2)gX'Z{M; 8*0(03), S(c)) 

= V^Z(M;0„S(c)), 

(2.23) 

where we have used the invariance of the topology of M under gluing a Tgn-i to 

M and the equations (2.9), (2.20) and (2.21). 

Next we will obtain the fusion rule of the states on S* . Define the following 
rank-three symmetric tensor: 

iV*^^^ = Z(5^;f,0^0'=), (2.24) 

where the manifold truly treated can be uniquely obtained by taking away 3 D- 
simplices from S . From now on, we will use such implications for simplicity. Using 
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this tensor (2.24), one can reduce the number of S s in correlation functions: 

Z{M-i:{c),(^')Ni^^ = Z{M-i:{c),BQ{c^)){BQ)\ci){BQUc2)g'i^ 

= Z{M- S(c), i?o(ci))^^7Z(5^; i?o*(c2), 0^ 0^) 
= Z{M] S(c) 'i>' ^^^ 



M y 5 



(2.25) 
Here, from the third to the last line, we glued the manifold obtained by taking 
away one D-simplex from S and the manifold obtained by taking away one D- 
simplex from M. Since the former manifold is topologically one D-simplex, the 
resultant manifold is simply M. In this derivation, we have used the constraint 
(2.9) for the gluing operation, (2.17), (2.18), (2.21) and (2.24). Taking particularly 
the three-sphere function as the partition function in the left-hand side of (2.25), 
one obtains 

Ni^^Nk^"^ = Ni^^"^ = Ni^^Nj"^, (2.26) 

where A^j-' "^ is the four-sphere function defined as same as (2.24). 
By the similar discussions as to derive (2.25), we obtain a formula 

Z{M- S(c), (t>')Z{S^; (t>^) = Z{M; S(c)). (2.27) 

This formula implies also that, if the TLFT is not null, the dimension of the 
physical Hilbert space on S must be non-zero. Substituting the first partition 
function of the left-hand side with a three-sphere function, we obtain 

2.28) 

where we have used the fact that S with two holes is topologically equivalent to 
TsD-i and (2.19). 
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From (2.25), (2.26) and (2.28), one can see that the physical states 0*s form 
a commutative algebra with an identity. This algebra might be decomposed into 
mutually annihilating sub-algebras, and the decomposition will give the decompo- 
sition of the TLFT. One can discuss the decomposition of the algebra in general, 
but in the present paper we restrict our discussions only to the case that the algebra 
has no radical . This is because, if a radical exists, one will have to discuss a wider 
class of TLFTs than those with (2.4). Let (pR be a radical element, and consider 
a volume-dependent TLFT Zs^^M; A; S(c)) = Z{M; S(c), exp{A6r(f)ji)). Since (pR 
is a radical element and a physical state, the Zg^ is a polynomial function of A. 
Thus, to discuss the neighborhood of a TLFT with radical, it is not sufficient to 
consider only the volume-dependence of (2.4). 

Since an algebra without radical can be decomposed into a direct sum of simple 
sub-algebras, the commutative algebra can be simultaneously diagonalized. We call 
the following the canonical basis ' ''' : 

r^y = ^y, r]ij = 6ij, m^ = \'6'^6'^ (A* ^ 0), (2.29) 

which is supposed to be used from now on. From (2.25) and (2.29), we obtain the 
fusion rule: 

0V^' ~ 5'^\'(t)\ (2.30) 

One can see in (2.28) that the one-sphere function is the inverse of the coeffi- 
cient of the three-sphere function: 

Z(5^;0*) = l. (2.31) 

From (2.27) and (2.31), we are lead to the fact that any correlation(partition) 



* We are not sure whether a TLFT can have an algebra with a radicaL In two dimensions, 
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no such TLFTs are known ' . One can reaUze such cases only by taking a certain infinite 



limit 
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function can be expressed in a sum 

^ 1 
Z{M-Mc)) = J2-Z\M-Mc)), ^^^^^^ 

where A^ is the dimension of the physical Hilbert space on S . In other words, 
(2.32) is 

^ 1 
E Y^^' - 1- (2-33) 

From now on, to avoid confusions, any repeated indices of physical states without 
the summation symbol are not supposed to be summed over. 

To show that (2.32) is a direct sum, we will prove the orthogonality of Z*s. 
Suppose M3 is obtained by gluing Mi and M2. Then, by using (2.9), (2.30) and 
(2.32), 

J]lz^(Mi;Si(ci))(?,^"^-i-ZJ-(M2;S2(c2)) = -i-Z(M3; £3(03), ^^ ^^ 



C-inn 



'; V ; 



= <5*^'i-Z^(M3;S3(c3)). 

(2.34) 

We will next discuss the decompositions of the physical Hilbert space and the 
physical correlation(partition) functions. The projection operator specifying the 
physical Hilbert space can be decomposed into 



Pse = $^P|:e 



PI,^P^^ = 6^^P^^ (2.35) 

^I:e(c2,ci) ^ ^gf'Z^i: X [0,l];S*(4),S(ci)), 

where we have used (2.12), (2.32) and the orthogonality (2.34). Define the i-th 
physical Hilbert space 7i|, ^*'* as the invariant subspace of the projection operator 

16 



iphys _ Y^ ^q_,phys,i 



Py,y.- Then (2.35) implies that H^ = Xlj (BH^ ^ ' . In particular, the dimension 
of each physical Hilbert space on S is one: dim(7Y^J'-f ) = 1- 

The physical states in the physical Hilbert spaces labeled by different i and j, 
Ti,^ ^^'* and Ti^,^^'"' [i ^ j), can not interact with each other. This is because the 
labeled partition(correlation) function is orthogonal to the physical Hilbert space 
with a different label: 

Z^(M;Si(ci),(S2y) = Z^(M;Si(ci),S2(c2))P4s.(c2,C3)(S2y(c3) ^ ^ 

' ' (2.36) 

= 5*^Z^(M;Si(ci),(S2)^'), 

where (S2)-' G Ti-j^ '^^''' , and, in this derivation, we have used (2.17), (2.35) and the 
orthogonality (2.34). The discussion is similar for the case of the dual states. Thus 
the linear map defined with the partition function Z in (2.16) is a direct sum of 
the linear maps defined with jtZ^s, which are the maps from T-Cj^ ^^'* ® ■ ■ ■®Ti-^ ^^'^ 
to 7Y|; ^^'* ■ ■ ■ H-^ ^^'*, respectively. This concludes that the original TLFT is 
a direct sum of the theories whose partition functions are given by j^Z^s. 

The orthogonality (2.34) implies that jtZ'^s of each i satisfy the defining con- 
straints of a TLFT (2.8) and (2.9). Thus jjZ^M; S(c))s of each i define a TLFT. 
Since the physical Hilbert space on S* of each TLFT is one-dimensional, each 
TLFT is an irreducible TLFT, which cannot be decomposed anymore into a direct 
sum of TLFTs. 

In the present paper we do not discuss the general case of the algebra, but one 
will be able to extend the discussions so far to the general case. An irreducible 
TLFT will have the algebra which cannot be decomposed anymore into a direct 
sum of mutually annihilating sub-algebras. 

Decomposing the TLFTs (without radical) in (2.10) into a direct sum of irre- 
ducible TLFTs, a partition function of a volume-dependent TLFT is a direct sum 
of those of irreducible TLFTs with volume-dependent numerical factors: 

N 
Z{M; A; S(c)) = J] Zi{M- S(c))(cu,)^, (2.37) 

i=l 
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where the Zi has only one physical state on S ^ . 

Since a volume-dependent TLFT with (2.4) is a direct sum of irreducible 
TLFTs, the natural definition of a physical Hilbert space of a volume-dependent 
TLFT is the direct sum of the physical Hilbert spaces of the irreducible TLFTs: 
^p ys _ ^^^^ ©7i|, ^^'*. In fact, a partition function of a volume-dependent TLFT 
can be regarded as a multi-linear map among Ji^ ^* as same as the map (2.16) of a 
TLFT. The physical correlation function defined similarly to that of TLFTs is in- 
dependent of the simplicial decomposition of M and its boundary, but is dependent 
on the volume of the manifold M. 

2. 5 Renormalization group flow 

Since a physical correlation function of a volume-dependent TLFT is dependent 
on the volume, there is a renormalization group fiow. A renormalization group 
flow of a lattice theory is determined by how the change of the lattice structure 
is absorbed in the change of the parameters of the theory without changing the 
physical outcomes. Since a physical correlation function of a volume-dependent 
TLFT depends only on the volume A of the lattice M, the invariance of the physical 
outcomes imposes only that {uJi) = (uJi) for each i under the change of the 
volume from A to A' . Taking the extensive variables pi = ln{uji) as the parameters 
of the theory, the flow is piA = p'^A' . This implies the volume-dimension of the 
parameter pj is 1, and this agrees with the Nienhuis-Nauenberg criterion 

The flxed points of the renormalization group flow are the points with Ui = 
0, l(or Pi = — oo, 0) for all i. The point pi = — oo is the infrared flxed point, and the 
point Pi = is the ultraviolet flxed point. Since the volume-dependence vanishes 
on such a flxed point, the theory becomes a TLFT. On the other hand, one can 
make a volume-dependent TLFT from a TLFT by (2.37) with certain cj^s. The 
original TLFT corresponds to the point uji = 1 for all i. Hence a TLFT can be 
regarded as being on a flxed point in the parameter space of a volume-dependent 
TLFT, and this flxed point has relevant operators with volume-dimension one. The 
number of such operators is equal to the dimension of the physical Hilbert space 
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on S^-'^ of the TLFT. 

2. 6 First order phase transitions and order parameters 

We will consider the thermo-dynamical limit A -^ oo. From (2.37), a partition 
function of a volume-dependent TLFT is 

iV 

Z{M- A) = Y, UM) exp{p,A), (2.38) 

where M is without any boundaries and A^ is the number of the irreducible TLFTs. 
The free energy per /^-simplex in the thermo-dynamical limit is 

1 
/ = — lim — InZ = — maxj(pj), (2.39) 

where maxj(pi) is for taking the pi which has the largest real part among all the 
PiS. Thus the point Pi = for all i is a multiple first order phase transition point 
with N different phases around it . In general, a TLFT is on a multiple first order 
phase transition point, and the number of the phases around it is equal to the 
number of the irreducible TLFTs of which the TLFT is the direct sum. 

The order parameters are appropriately given by using the one-sphere func- 
tions. Define the one-sphere expectation values as follows: 

where 0* is a physical state on S* , and we take the canonical basis to label 
them. Using the orthogonality (2.36), we obtain 

where we assumed the system is in the m-th phase, that is, pm has the greatest real 
part among pi{i = 1, ■ ■ ■ , A^). Thus the one-sphere expectation values give good 
order parameters. 



• In general, Zi{M) might be zero, and then the discussions must be changed. Here we set 
aside such special cases. 
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Comments are in order. So far we have assumed implicitly that there are no 
constraints on the parameters of the theories and that we can take any values for 
the parameters cjjs or piS. But in some physical cases we cannot. An example 
is the 2-dimensional lattice QCD in the following subsection. This theory is a 
volume (area)-dependent TLFT, but cannot have any first order phase transition 
points if the positivity of the local weight is assumed. 

2. 7 An example — 2D lattice QCD 

As a simple but physically interesting example of a volume(area)-dependent 
TLFT, we consider the 2-dimensional lattice QCD ' ' ' "" . As the index set Xi 
for l-simplex(edge) we take the Lie group G = SU{n),U{n). The index set is 
infinite for this model, but nonetheless we can treat it straightforwardly in the 
same way as in the preceding subsections. The index set for O-simplex(vertex) is 
not considered in this model. The partition function for a single 2-simplex(triangle) 
and the gluing tensor are defined as 



Z{T^; l;dT^{gi,g2,g3)) = Cg,,g^,g, = z{gig2gz), 

z{ugu-^)=z{g), (2.42) 

9'''''=5igi,g2'), 

where u,g,gi G G, and T^ is a 2-simplex(triangle). We take the local weight z(g) 
as a positive real function of the group element of G, and it satisfies the invariance 
in the second line. The g9^92 jg ^]-^g gluing tensor for a pair of 1-simplices, and is 
expressed by the delta function associated to the Haar measure of G, which is the 
measure used in summing over the internal colorings. 

To show that the elementary data (2.42) define the partition functions of a 
volume-dependent TLFT, the following invariance of the product of two local 
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weights under the flip move (Fig. 4) is essentiah 



dgz{gig2g)z{g ^g^gA) = / d^g^ W)^(^i^2(^2 ^9gA))z{{g2 ^ggd ^gsgA) 

dgz{giggA)z{g^^g2gz)- 



(2.43) 
It is known that all triangulations of a two-dimensional surface with the same 
number of triangles can be generated by the flip move . A partition function 
Z[M; A; S(c)) can be obtained from the elementary data (2.42) by using the rules 
(2.2) and (2.3). The invariance (2.43) guarantees that the obtained partition func- 
tion is in fact independent of the triangulation of M and that the partition functions 
satisfy (2.2) and (2.3). 

To perform the Laplace decomposition of the area-dependent TLFT, we will 
calculate Z(B'^; A;dT'^{gi, g2, gs)), where B^ is a two-dimensional ball. The in- 
variance in the second line of (2.42) enables one to expand the local weight z in 
the characters of G: 

z{g) = J2^RdRXRig), (2.44) 

R 

where d^ is the dimension of the irreducible representation R of G, and A^^ is the 
coefficient of the expansion. The partition function Z{B'^; A;dT'^{gi,g2,g3)) can 
be calculated easily using this expanded form (2.44) and the orthogonality of the 
characters: 

Z{B^;A;dr\gi,g2,g3)) = 5](Ai?)^rfi?Xi?(^i^2^3). (2.45) 

R 

This result (2.45) and the orthogonality of the characters imply that each d^XR 
gives the local weight of a TLFT. To see that each weight defines an irreducible 
TLFT, we will check dim(?i^i^* ) = \. We take a three-gon as the simplicially 
decomposed S^. The projection operator from a three-gon to another is calculated 
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to be 

= XR{{9l929^)-')XR{9'i9'29'i), ^^.46) 



This projection operator projects any wave-function 4>{gi, g2, 93) to XR{{9l929i) ^)- 
Thus 'H^gi'^ = {XR{{9l929i)^^)} and its dimension is one. 

Since the 2D lattice QCD is an area-dependent TLFT with the decomposition 
as a direct sum of an infinite number of TLFTs, one might expect that it has first 
order phase transition points. But, in fact, it is not vahd because the coefficient 
Aq for the trivial representation in (2.44) has always the greatest absolute value. 
This inequality can be shown as the following. Since the local weight z is assumed 
to be positive real, 

|Ai?| = 1^ / dgxR{9)z{9)\ = 1^ / d9^Pi{9)z{9)\ < / dgz{g) = Aq, (2.47) 
where piS are the eigenvalues in the representation R. 

3. Neighborhood of topological lattice field theories 

In the previous section we investigated volume-dependent TLFTs and showed 
that TLFTs are on fixed points of the renormalization group fiow as well as on mul- 
tiple first order phase transition points in the parameter space of volume-dependent 
TLFTs. In this section, we will consider general perturbative deformations of 
TLFTs. To treat such deformed systems, we will introduce a perturbative method, 
which is a kind of cluster expansion . We assume that such a perturbative method 
is available in the neighborhood of a TLFT. 

Consider a manifold M simplicially decomposed with a certain number of D- 
simplices. The partition function of M of a TLFT is expressed as a product of 
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local partition functions of the TLFT: 

C q b (3.1) 

Z(r^(g);c) = Z(r^(g);9T^(c)), 

where the summation is over the inner colorings, and the products are of the local 
partition functions over all the D-simplices and the gluing tensors. 

We consider general changes of the local partition functions; Z{T {q)]c) -^ 
Z{T^{q)]c) + 5Z{T^{q);c). We assume that 6Zs are small enough to justify 
the perturbative treatment in the following subsection. We will not consider the 
changes in the gluing tensors for simplicity, because they will not change the dis- 
cussions essentially. Since we consider general changes, the deformed system is not 
a TLFT anymore. Thus the partition function depends on the simplicial decom- 
position of the manifold. Thus, from now on, we are supposed to consider certain 
simplicial decompositions and will take the thermo-dynamical limit that the total 
number of the D-simplices in M is taken to infinity. 

3.1 Cluster expansion 

Here we will introduce a systematic perturbative method, which approximates 
the partition function and the correlation functions of the deformed system. This 
method is a kind of cluster expansion , but is unusual in the sense that it is a 
cluster expansion of an operator. 

Here we begin with the partition function. The partition function of the de- 
formed system is defined by 

Z{M; SZ) = J2 ll(^ir''{qy,c) + 8Z{r^{q)- c))\{gl^\ (3.2) 

c q b 

where q is the label for the D-simplices in M. Expanding this partition function 
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(3.2) in 6Z, we obtain 



A 



n=0 ' qi,--,q„ ^ ' ' 

UQ = Z{M), 

where A is the total number of the D-simphces in M, and Z{M) is the partition 
function of the original TLFT. The Un{qi,- ' ' ^Qn) describes an 0{{6Zy'') contri- 
bution to the partition function, and the g^s specify the n 6Zs coming from the 
D-simphces qi,- ■ ■ ,qn- The Un{qi, ■ ■ ■ , q-n) is supposed to be non-zero only when 
gi, ■ ■ ■ , gn are different from each other. 

If we assume that the simplicially decomposed manifold M has an enough 
'extension', the lower order contributions can be calculated using the correlation 
functions of the original TLFT on M: 

Mn(gi, ■■■,qn)= Z{M; Uniqi, ■ ■ ■ , qn)) + 0{n - n{M))Z{Mi; U'^{qi, ■ ■ ■ , g^)) + ■ ■ ■ , 

(3.4) 
where n{M) is the number of D-simplices below which one does not have to con- 
sider correlation functions on manifolds different from M, and 6{n) is the discrete 
analogue of ^-function: 1 for n > and for n < 0. An example of M and qi,- ■ ■ ,qn 
that one has to consider correlation functions on manifolds different from M are 
given in Fig.5. Uq is the identity and Un takes non-zero values only when the g^s 
are different from each other. Since n{M) = 1 in D = 1, we will consider only 
the D > 1 cases. The D = 1 case can be treated very easily in another way. We 
assume that the thermo-dynamical limit is taken by refining the lattice structure 
of M almost uniformly. Then n{M) ^ oo in the thermo-dynamical limit, so that 
the right-hand side of (3.4) will be well approximated by the first term. Thus we 
will consider only the first term from now on. 

The definition of the first term of (3.4) is as follows. Consider a lattice M' 
obtained by taking away n D-simplices gi, ■ ' ' 5 9n from M. Then M' is a lattice 

24 



with some boundaries S, and a partition function Z{M'; S(c)) of the original TLFT 
can be associated to the lattice M' . Contracting the colorings of Z{M'; S(c)) with 
those of 5Z{T (qi); c)s, one obtains Z{M; Un{qi, ■ ■ ■ , qn))- Thus Un{qi, ■ ■ ■ , q-n) 
represents the operation to take away n D-simplices qi,- ' ' ^Qn from M and contract 
the colorings of the boundaries with those of 5Z{T (qi); c)s. 

If the n D-simplices qi^- • • ,qn are separated, Z(M; f/„(gi, ■ ■ ■ , g„)) is given by 
contracting the colorings of the n-sphere function of the original TLFT with those 
of 5Z(T {qi);c)s. Since the n-sphere function is independent of the positions qi, 
the qi dependences of Z{M; Un{qi, ■ ■ ■ , qn)) come only from those oi5Z{T^{qi); c)s. 
This property simplifies the systematic calculation of (3.3) drastically, because the 
leading contributions come from the cases that the gjS are separated. 

To consider sub-leading contributions, we introduce new operations Wn{qi, ■ ■ ■ , qn) 
satisfying 

(s fci times kn times 

EikA ^n 
I— 1 * 

(3.5) 
where an denotes the permutation maps {l,---,n} -^ {l,---,n}. The relation 
(3.5) defines Wn iteratively with Un and Wi ~ Wn-i- The relation between WnS 
and Z{M; 6Z) is given by 



Z{M;6Z) = Z{M;exp{W)), 

n=l qi,---,qn 



(3.6) 

nj, 



which can be proven formally using the definitions (3.3), (3.4) and (3.5). 

The relation (3.5) is of course meaningless before the products of WnS are 
defined. From (3.5) one obtains Wi{q) = Ui{q) without any such products. But 
to obtain W2 we have to define the product Wi{qi)Wi{q2) ■ In fact there are 
various ways to define the product. If gi and g2 are separated, a natural definition 
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of the product would be to identify Wi{qi)Wi{q2) with U2{qi,q2), because the 
former can be regarded as the operation of taking away the Z)-simphces gi and 
g2- But how about the cases that gi and g2 coincide or are next to each other? 
We use 'operator sphtting regularization' for the definitions of the products in 
general. The meaning is as follows. For example, Z{M;Wi{qi)Wi{q2)) is defined 
by contracting the colorings of 5Z(T (qi);Ci) (i = 1,2) with those of the two- 
sphere function Z {M ; dT (ci) , dT (02)) of the original TLFT. This definition 
coincides with U2{qi, ^2) if qi and g2 are separated. From now on, we assume that 
such 'operator splitting regularization' is used as the definitions of the products of 

WnS. 

Since W2{qi,q2) = U2{qi,q2) - Wi{qi)Wi{q2) from (3.5), 1^2(91, ^2) takes non- 
zero values only when qi and g2 coincide or next to each other. Such a cluster 
property of WnS can be proven in general as follows. Consider two sets of positions: 
Q = {qi, ■ ■ ■ , qm}, Q' = {o'm+i, ■ ■ ■ , Qm+m'}- Suppose these two sets are separated 
from each other. We will prove by iteration that Wm+m'{qi,- ■ ■ ,qm+m') is zero 
always for such separated Q and Q' . We start with the following equation which 
is valid because the two sets are separated: 

Um+m'iqi, ■ ■ ■ , qm+ni') = Umiqi, " " " , qmWm'iqm+l, " " " , qm+m'), (3.7) 

where the definition of the product of f/mS follows that of WnS. As was shown 
above, the cluster property is valid for m = m' = 1. Assume that m + m' > 2, and 
that Wn {n < m + m') satisfies the cluster property. Then, from the relation (3.5), 

Um+m'{qi, ■ ■ ■ 5 qm+m') — Wm+m'{qi, " " " , qm+m') 

(m+m' — I ^ \ 

n ki\{i})k^ I 5Z ^l(^a„+„'(l))---^'n+m'-l(---,?a„+„,(m+m')) 

Em+ra' —1 . 

z — 1 * 

= Umiqi, ■ ■ ■ , qm)Um'iqm+l, " " " , Qm+m'), 

(3.8) 
where we used the cluster property of Wn {n < m + m,') from the second to the 
last line. Comparing with (3.7), we obtain Wm+m'ill, ■ ■ ■ , Qm+m') = 0. 
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Comments are in order. Because of the topological nature of the partition(correlation) 
functions of the TLFT, in the 'operator splitting regularization', one has to take 
care only of the topology of the neighborhoods of the operators. Thus it is rather 
simple in lower orders, but is complicated in higher orders since the topology of 
the neighborhood becomes complicated. But we are allowed to choose the neigh- 
borhood of an operator to make calculations as simple as possible. The difference 
of the choices will be absorbed in the operators with orders higher than the con- 
sidering order. Consider, for example, a product Wn{qi, ■ ■ ■ , (ln)Wi{q), where the 
D-simplices qi,- ■ ■ ,qn surround an area of the lattice. Naively there are two possi- 
bilities of the splitting, say the inside or the outside of the area. We should define 
the inside of the area as a part of the neighborhood of gi, ■ ■ ■ , g„, because then the 
splitting is always done to the outside and this facilitates the calculations. 

The partition function Z{M;J]{c);6Z) with colored boundaries S(c) can be 
also calculated by the cluster expansion. We start with the analogue of (3.3) and 
(3.4): 

Z{M; S(c); 6Z) = Z{M; f/^^")), 

n^Oqi,---,qn 



U ^^' = y > ^t/n-igi, ■■■,?«;, (3.9) 



U^^^^ = S(c). 



Special cares must be taken if some of the qi,- ■ ■ ,qn are next to S. Define W^ to 
single out such boundary effects: 

n ^ 

Un " (gi, ■■■,qn)= ^Yl rr,\(r, - rr.W ^'^ ^ (^'^"(1)' ' ' ' ' ?a„(m))f^n-m(ga„(m+l) , " " " , Qa^in)), 

1 1 L t\ I L libit 

m=0 o-„ ^ ' 

(3.10) 
where Un{qi, ■ ■ ■ , qn) is defined by splitting operators from S. One can show that 
Wn satisfies a cluster property, that is, Wn {qi, ■ ■ ■ , qn) takes zero if there 
exists a non-empty subset Q of {gi, ■ ■ ■ , qn} separated from its compliment Q as 
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well as from S (Fig. 6). Using the Wn , the partition function is expressed as 
Z{M; S(c); 6Z) = Z{M; W^^^'^ exp(W^)), 

^"^^^ = E E >"^^^fe,---,^n), (3.11) 

n=Oqi,---,qn 

<(^) = S(c), 

where the definition of W is given by (3.5) and (3.6). 

Consider the case that S is composed of two boundaries Ej (i = 1, 2) and 
they are separated by d lattices. Consider a set Q = {qi,- ' ' yin} for which 
^n{<li,- ' ' jQn) 7^ 0. Then, if n is smaller than the distance d, the set Q is a 
sum of two sets Qi and Q2, where Qi is in the neighborhood of Si but is sepa- 
rated from Q2 and the boundary S2, and vice versa for Q2 (Fig. 7). One can prove 
that, if n < d, Wn can be factorized: 

n 

= J2Y1 m\(n-m)\ ^"'' ""' (^<^4l)' ' ' ' ' 9'T„(m))W"n-m (?a„(m+l), " " " , ga„(n)), 
m=0 cr„ '^ ''■ 

(3.12) 
where W^^{qi,- ■ ■ ,qn) {i = 1,2) has the cluster property that it is zero if there 
exists a subset of {gi, ■ ■ ■ , g^} separated from the compliment and the boundary Sj. 
An important point is that W^"- is independent of the other boundary. Substituting 
(3.11) with (3.12), we obtain 

Z{M; Ei2(c); (5Z) = Z{M; (i^Si(cO|^S2(c2) ^ w^^'^)^ exp(l^)), (3.13) 

where W^ ^^ is the order of {5Z) . 

3.2 First order phase transitions 

In the previous section, we investigated volume-dependent TLFTs, and showed 
that a TLFT is on a multiple first order phase transition point and that the number 
of the phases around it is equal to the number of the irreducible TLFTs of which 
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the TLFT is the direct sum. This statement is exact because no approximation 
is used, but is vahd only in the parameter space of a volume-dependent TLFT. In 
this section, we will show, in the approximation of the first order of the cluster 
expansion, the same statement in the parameter space of the general deformations 
of the local weights. 

The operator Wi{q) is to take away a D-simplex and contract the colorings 
with those of 5Z(T (q);c). Thus the 'operator splitting regularization' implies 

n n 

Z{M- S(c), n Wiiqi)) = Z{M; S(c), i?o(ci), ■ ■ ■ , i?o(cn)) J] ^^'^^(^''te); cD 

n 

= Z{M- S(c), 0^S ■ ■ ■ , 0^") n <PjAcr)9t^Z{r^{q,)- cO, 

(3.14) 
where Bq = dT and 0iS are the physical states on S* . Here we have used 
(2.21). Thus, in the first order of the cluster expansion, W is 

W = ^Wi{q) = A6Zi(f)\ 

1 (3.15) 

SZi^jJ2<l'^ic)gi6ZiT^{qy,c'), 

g 
where A is the total number of the D-simplices in M, and SZiS are the mean values 
of the deformations of the local weights projected to the physical states on S . 

Using (3.6), (2.30), (2.33) and (2.32), Z{M; 5Z) can be calculated easily in the 
canonical basis: 

Z{M-5Z) = Z(M;exp(W^)) 

= J]-Z(M;0*exp(W^)) 



1=1 

N 



j=l n=0 

J]-Z^(M)exp(AAMZ,), 



(3.16) 



.=1^ 
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where N is the dimension of the physical Hilbert space on S , or the number 
of the irreducible TLFTs of which the TLFT is a direct sum. This is formally 
equivalent to (2.38) with pi = X^6Zi, and again the TLFT{6Zi = 0) is on a multiple 
first order phase transition point with N different phases around it. 

In the previous section, the order parameter is given by the one-sphere expec- 
tation value. Analogously, here we consider the following expectation value: 

Using the discussions of the cluster expansion with boundaries in the previous 
subsection, we obtain, in the same way as (3.16), 

iV 
Z(M; 6Z; Bo{c))^\c) = J^ Z{M; (0^ + O(5Z)^0^) exp(W^)) 

';' (3-18) 

= X^('^^ + 0{5Zfj)Z{M; 0^') exp{A\^Zj), 
i=i 

where 0{5Z) ■ comes from the correction to the boundary Bq as in (3.11). Thus, 
substituting (3.17) with (3.16) and (3.18), we obtain 

>^.k.M,5Z ^ li^ 1 Ef=i(^'i + 0{6Zfj)Z{M- <pn expjAX^Zj) 
A^ ^"^A^ J2f^^ijZ{M-<pj)exp{AXJSZj) (3.19) 

= 5>^m + Oi5Zf^, 

where we assumed that the system is in the m-th phase, that is, W^SZm has the 
greatest real part among ySZi {i = 1, ■ ■ ■ , N). 

3.3 Correlation lengths of fluctuations associated to trivial local operators 

We here discuss the correlation lengths of fluctuations associated to trivial 
local operators in the cluster expansion and show that they vanish in the limit of 
TLFT 6Z -^ 0. Consider a closed lattice M and make a new lattice M' with two 
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boundaries Si^2 by taking away some of the D-simplices in M . We assume that 
one of the boundaries Si is topologically trivial, that is, Si can be surrounded by 
a S in M' and separated from S2 by the S . Thus what we consider here 
are the fluctuations associated to such a topologically trivial local boundary Si. 

Let Oi,2 be the operations to make the colored boundaries Si^2(ci,2), respec- 
tively. As was discussed previously, these operators are modified in the cluster 
expansion. Assume that each operator is modified to O'^ 2 respectively. Then the 
composite operator O1O2 is modified to 0'iO'2 + O^^ where the O^ is of order {6Z) 
when the two boundaries are separated by d lattices. Then the expectation value 
of the correlation function of the two operators is calculated as 

u^ .^ .M,5Z ^ y Z{M-0^,02;5Z) ^ f ZiM^O[0^eMW))_ ZiM;OsexpiW)) \ 

^uiU2) -^^^ Z{M;6Z) A^oc\ Z{M;exp{W)) Z{M;exp{W)) J 

(3.20) 

where W is that in (3.6). Since the operator Oi is a trivial operator, the first term 

can be factorized: 

j.^ Z{M;0[0',exp{W)) ^ ^.^ Ek=l Z{S^; <l>kO'i)Z{M; O'^ct^^ exp(iy)) 

^hmZ(^^;0^O;)?^^^^^^-P(^)), 
A-00 ^ '^"^ ^^^Z(M;exp(W^)0-)' 

(3.21) 
where we cut M' at the S* and used the constraint (2.9), (2.21) with S = S* 
and (2.33) in the first line, and we assumed that the system is in the m-th phase. On 
the other hand, the expectation values of the one-point function of each operator 
are 

(Oi)^'^^= hm Z{M-exp{W)r)Z{S^-A'^0[) 



A-00 -^Z{M- exp{W)<P^) 

= lim X^Z{S^-(rO[), (3.22) 

A^oo 

ip2)^^'^= lim ^(^;^2exp(W^)0'-) 



A^oo Z(M;0"^exp(iy)) ' 
where we have used the same cut as above in the first line. Thus, combining (3.20), 
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(3.21) and (3.22), we obtain the connected correlation function as 

A~*oo z(M;exp(M/)) 

This concludes that the correlation length is of order — l/ln((5Z), which vanishes 
in the hmit of TLFT (5Z ^ 0. 

3.4 An example — 2D Q-state Potts model 

The previous example of 2D lattice QCD cannot be on a multiple first order 
phase transition point, if we respect the physical constraints. But the 2D Q-state 
Potts model has a multiple first order phase transition point of a TLFT in the 
physical parameter space 

The Potts spins are located on the 1-simplices, and they interact among nearest 
neighbors and with external magnetic fields. This model is defined by the following 
local partition function and gluing tensor: 

g'^ = 6'^, djk = KW{t,j)W{j,k)W{k,t), 

(3.24) 
W{i,j) = exp{p{6ij -l) + hi + hj), 

where [3 is the inverse temperature, hiS the external magnetic fields, and the Roman 
indices run from 1 to Q. As was discussed in the previous example, this model will 
be a volume-dependent TLFT if the 'flip invariance' is satisfied. One can easily 
show that the necessary and sufficient condition of the invariance is /3 = or 00 

For /5 = 00, W{i,j) = 6ij exp{2hi). To perform the Laplace decomposition, we 
calculate Z(B'^; A; _Bo(«, j, k)), where B^ denotes a disk and Bo(i,j, k) a three-gon 
with coloring i, j, k: 

Q 
Z{B''; A; Bo{i,j, k)) = 6ij6ik («:exp(6/i,))^ = Yl ^'^^f^k {'^exp{6ha)f . (3.25) 

Q=l 

The local weight 5f 5^5^ trivially defines an irreducible TLFT for each a. Thus 
the point K,exp[6ha) = 1 for all a must be a multiple first order phase transition 
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point with Q different phases around it. On the other hand, since, after redefining 
K and the magnetic fields properly, this point corresponds to the point of zero 
temperature and zero magnetic fields of the Q-state Potts model, it is a multiple 
first order phase transition point with Q different phases labeled by the directions 
of the spin. Thus the two results coincide. The cluster expansion of the partition 
function near the TLFT has a direct relation to the low temperature expansion of 
the Potts model. This suggests that the cluster expansion near a TLFT is at best 
an asymptotic expansion in general. 

For /9 = 0, W{i^j) = exp(/ij + hj). One obtains easily 

Z(S2; A; Bo{t,j, k)) = «;(Eji^ eM^hi)f/Y{^tl exp(4/iO)~^/' exp(2(/ii+/i,+/ifc)). 

(3.26) 
Here one can see that the local weight (X]l=i exp(4/;,;))~^/^ exp(2(hi + hj + hf^)) is of 
the form i^iLpjLpk with Lpi = (^Jl^ exp(4/i;))^ ' exp(2/ij). Hence the local weight 
can be transformed to the form 6l5^5j^, and it defines an irreducible TLFT. This 
concludes that the TLFT at /? = and k{Y^^^^ exp(4/i/))^'^ = 1 is an irreducible 
TLFT and is not on a first order phase transition point. On the other hand, this 
point corresponds to the infinite temperature limit of the Potts model, and hence 
is not on a first order phase transition point. The cluster expansion near the TLFT 
corresponds to the high temperature expansion of the Potts model. 



4. Neighborhood of three-dimensional Zp model 

So far we have discussed the physical states on S of a TLFT and their 
relations to the phase structure of first order phase transitions. We gave some 
examples in two dimensions. Since, in two dimensions, a boundary has always the 
topology of S^, all the physical states are those on S^. Thus only the structure 
of first order phase transitions is associated to a two-dimensional TLFT. But how 
about the physical states on non-trivial boundaries in higher dimensions? We will 
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not give here any general statements about it. We will discuss it in a simple three- 
dimensional model on a manifold with a certain topology. The method is again 
the cluster expansion. But here higher orders will play essential roles. 

4-1 Definition and physical states 

The model we discuss here is a simpler analogue of the model by Turaev and 

foil 

Viro . We consider the Zp gauge group in place of the quantum groups. (See [22] 
and [25].) 

Consider a simplicially decomposed three-dimensional manifold M possibly 
with some boundaries. The degrees of freedom are on its 1-simplices, and they 
are integers running from to p — 1. The coloring of the lattice is performed by 
the assignment of these integers to all the 1-simplices. The weight for each inner 
0-simplex is 1/p, and that for each 0-simplex on boundaries is 1/^^- The weight 
for each 3-simplex is given by the Zp analogue of 6j-symbol: 



^' -^2 JS (p) (p) (p) (p) 

"7i+i9+r^.0"fi+iV,-ifi.0"i9-i4+'/fi.0"is+74-7V,.0' y^--^) 

J4 J5 J6 



dp) 



where 6^^ takes 1 for i = j mod p, but otherwise. The weight (4.1) is a flatness 
condition for the gauge group Zp, that is, the summation of jjS over the 1-simplices 
surrounding each 2-simplex should vanish with mod p. The partition function is 
defined as 



Z{M; E) = p-(«>+«''/2) y^ Y\ 



j T 



■T -T -T 

h h h 

■T -T -T 

M h h 



(4.2) 



where Vi and Vh are the total numbers of the inner 0-simplices and those on the 
boundaries, respectively. The summation is over all the colorings, that is, all the 
assignments of the p integers to the 1-simplices. jj represents the six integers 
associated to each 3-simplex and the product is taken over all the 3-simplices in 

M. 

All the triangulations of the inside of a three-dimensional manifold can be 
related by the two local moves (1,4) and (2,3) and their inverse moves (Fig.8). 
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The partition functions (4.2) are invariant under these moves, and hence they 
define a TLFT. 

The boundary of a 3-manifold is a closed 2-nianifold. Let Bg denote a trian- 
gulated surface with genus g. Since (4.2) consists of the products of the fiatness 
conditions, one can show easily that the basis of the physical Hilbert space on 
Bg can be labeled by the currents on a and j3 cycles of Bg-. Tig = {{jajpi (i = 
1) ■ ■ ■ )fl')}|jai? J/3i = 0, ■ ■ ■ ,p — 1}. Hence the dimension of the physical Hilbert 
space is given by dim(7ig) = p^^ . Especially, the dimension of the physical Hilbert 
space on the sphere is one, so this system is not on a first order phase transition 
point. 

4-2 Evaluation of the partition function 

We will calculate by the cluster expansion the partition function of the system 
obtained by shifting the local weights by little amounts from those of the Zp model. 
Since the deformed system is not a TLFT anymore, we have to fix the lattice. The 
lattice we consider is a simplicially decomposed manifold Mg with the topology of 
S^ X Bg. Since we finally consider the thermo-dynamical limit, we can consider 
a surface composed of 2-simplices in Mg which cuts Mg almost perpendicularly 
to the direction along S^ at any place t(G S^). We assume the foUowings for the 
triangulated surface t x Bg (Fig. 9): 

(1) The holes are regularly ordered and the length of the circumference of each 
hole is about /i lattices. 

(2) The length of the cycle winding through two neighboring holes is about I2 lat- 
tices. 
(3) h < h. 

In the thermo-dynamical limit, the length of S^ and the extent of Bg will 
become so large that the partition function of the deformed system can be well 
approximated by a correlation function of the original TLFT on Mg-. Z{Mg\ 6Z) = 
Z{Mg;exp{W)). 

The first order of the cluster expansion of W is the operator taking away one 
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3-simplex. Thus the 'operator sphtting regularization' imphes 

W = AMo, (4.3) 

where 0o is the physical state on S'^, and /o is the order of 5Z. Let us consider 
the p-th order {p > 1). The p-th order of W comes from the p 3-simphces forming 
a cluster, li p < li, the p 3-simplices can be always surrounded by a S*^ in Mg. 
Taking such a S*^ of a small size, one can define the neighborhood of the cluster 
as the neighborhood of the ball whose boundary is the S*^. Then the 'operator 
splitting regularization' projects the operation of taking away the p 3-simplices to 
the physical state 0o- Thus the cluster expansions below the /i-th order can be 
included in the form (4.3). 

New operators arise when we consider the /i-th order. There appear the cases 
that /i 3-simplices form the non-trivial loop surrounding one of the g holes in Mg. 
Then the 3-simplices can not be surrounded by a S*^ in Mg. The simplest definition 
of the neighborhood of such 3-simplices is to choose a small torus which surrounds 
the 3-simplices around the hole. Then such operations are projected to the physical 
states on the torus by the 'operator splitting regularization'. Thus we obtain 

9 p-i 
W = A/000 + 2^ E E 4<^fc ' (4-4) 

k=l i,j=0 

where T is the number of the lattices along S^. The (p^^ denotes the physical states 
on the torus around the k-th hole in Mg, and is labeled by the two integers along 
the a and /3 cycles of the torus. The //^s are the mean values of the coefficients 
coming from the projections of the operations to 0^'': fj'j = ^ ^^ fiji^)- Thus they 
are the order of (SZ)^. Since, from (2.33), the first term of (4.4) is merely the all 
over factor of the partition function, we will consider only the second term. To 
calculate Z{Mg; exp(W^)), we begin with determining the 'matrix elements' of (p^^ . 
Consider a manifold Tg = [0, 1] x Bg with an insertion of the operator (f/^ . Then, 
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using similar discussions to obtain the physical Hilbert space before, the 'matrix 
elements' of ^^-^ are obtained from the flatness condition as 

1=1 

where ■?/'^ and iIjq ' are the physical states on the boundaries B* and Bg^ 

respectively. Here we took the normalization of the inserted operator 0^^^ such that 
there appear no numerical factors in the right-hand side of (4.5). The case that 
there are two insertions can be considered by taking the square of the matrix (4.5), 
and we obtain the fusion rule 

For simplicity, we will use tp^^ = (1/p) X]f=o exp(i27rj7/p)0^ in place of 0^'', where 
the i in i2Ti denotes the imaginary unit, and do not confuse it with the label for 
the states. Then the fusion rule and W become 

Vk Vk ^iV ^jj' Vk ' 

9 P-i (4.7) 



fc=l i,j=0 



where h^-s are the order of (SZ)^ . Then, using the fusion rule (4.7) and Z{Mg] Lp^-^-'^ , ■ ■ ■ , (^9^ 
= 1 obtained by taking the trace of (4.5), we obtain 

Z{Mg- exp(iy)) = Z(M,; \[ J^ vl ^MTh%)) = \{Y. ^MTh'i^)- (4-8) 

k—l i,j=0 k=l i,j—0 

Thus the free energy of the system is 

F{Mg-6Z) = -EME e^P(^4))- (4-9) 

fc=l i,j=0 
4-3 Phase structure 
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Since we assume 5Z <^ 1, if one takes /i -^ oo in the thermo-dynamical 
limit, the system will become uninteresting. Thus we fix /i and I2 in the thermo- 
dynamical limit. Since then the number of holes g is proportional to the extent of 
Bg^ it is natural to consider the free energy per hole and per lattice along S^: 



1 
T,5^"oo gT' 



f{6Z) = Jim -T^F{Mg] 6Z) = -maxij{hij), (4.10) 



where raa.Xij{hij) is to take the hij which has the greatest real part among all the 
hij (^5 J = 0, ■ ■ ■ ,p— 1), and we assumed for simplicity that h^- does not depend on 
the hole: h^, = hij. From (4.10), we conclude that the TLFT is on a /i-th order 
discrete phase transition point with p^ different phases around it. 

Comments are in order. Firstly, the number of the different phases around 
the TLFT would depend on the physical constraints of the system. Secondly, the 
order of the phase transition /i is not universal, since /i depends heavily on the 
local lattice structures around holes. Thirdly, if /i is large, the phase transition 
is very weak and hence the phase structure might become different because of 
possible large fluctuations. Finally, the order of the phase transitions is /i only on 
the point of the TLFT {6Z = 0) in general. An example is as follows. Suppose the 
free energy is given by —raa,x{x'^y^~'^ , x'^ y ) {q 7^ q')- Then a phase transition 
line is at a; = y. But the line is a flrst order phase transition line except the point 
X = y = 0, where the order is /i. 

4.4 Correlation length and order parameter 

In the preceding section, we showed that the correlation lengths of the fluctu- 
ations associated to a local trivial operator vanish in the limit of TLFT. In this 
section, we investigate the correlation lengths of the fluctuations associated to local 
non-trivial operators in the perturbed Zp model. We will show that the correlation 
lengths of the fluctuations associated to the operators which wind through two 
neighboring holes diverge in the limit of TLFT. 

We flrst discuss the correlation lengths of the fluctuations associated to the 
torus around a hole. Consider a triangulated torus Bi around the k-th hole and 
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an arbitrary boundary S obtained by taking away some 3-simplices from Mg. We 
assume they are separated by d lattices. We start with evaluating the following 
partition function in the cluster expansion of the /i-th order: 

p-l 
hm J2 Z{Mg- ^{c),Bi{c'y,6Z)^^{c') = lim J^ ^(^^5 (S(c)'D^-'' .v^f^' + Od) exp{W)) 

d i',j'=0 



hm Z{Mg; {^{cYD^j^^^^r + Od) exp(W^)) 
1 ^oo 






= ^lim^(D;-'^^Z(M,; S(c); 6Z) + Z(M,; Odexp(l^))). 

(4.11) 
Here the modifications of the operators by the cluster expansion are considered, and 
D]^ ; ■, are the coefficients obtained by projecting the modified -Bi(c) to the y^^j^s. 
Hence -0^"^^, •, = S^f, 6^^J +0{6Z). We assumed that, in the limit T ^ oo, /i^^,^ domi- 
nates among h^j {i,j = 0, ■ ■ ■ ,p-l), and used exp(W^) = 11^=1 '^^j=o V'k exp{Th'yj) 
and the fusion rule (4.7) in the last line. The operator O^ is the modification of 
the composite operator like W^ ^^ in (3.13), and is of order {5Z) . Substituting 
(4.11) with S = 0, we obtain another relation 



hm J2 Z{Mg; Biic'y, 6Z)^'l{c') = lim /^L,^(M,; 6Z), (4.12) 

1 — >0O ' ' 1 — >CX) 

c' 

From (4.11) and (4.12), we obtain 

(Mc)V^fc) =^lim^ Z(M,; 5Z) (4.13) 

= (S(c))^-'5^(v9^^')^-^^ + (Od)^--^^. 

Hence the correlation lengths of the fluctuations associated to Bi are the order of 
— l/ln5Z, which vanish in the limit of the TLFT {5Z -^ 0). 

From (4.12), one can see that the order parameters are given by the following 
one-93 expectation values: 

^^ij^M^.Z ^ jy,^ ^ ^^ ^ om- (4.14) 

Next we consider the fluctuations associated to the torus winding through 
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two neighboring holes. We denote the physical states associated to the torus by 
r^-'p where kl labels the two neighboring holes and ij the physical states (Fig. 10). 
Using the similar discussions as before, the 'matrix elements' of these operators 
are obtained from the flatness condition: 



m=l 



where we normalized the Fs such that there are no numerical factors in the right- 
hand side of (4.15). For convenience, we will use the operators 7^] = ^ F^^ exp{i27cqj/p) 
in place of F. 

One has to take into account the deformations of the operators in the cluster 
expansion. But, since the main contribution comes from the original operators, we 
will neglect the deformations. This simplification does not change the conclusions 
concerning the correlation lengths. 

By using the fusion rule (4.7) and taking the trace of the product of the matrices 
(4.5) and (4.15), the one-7 function is evaluated as 

Z{Mg; 7i]; 5Z) = Z{Mg; 7^] exp(l^)) = J^ Z{TB,;^pt''-^\ ^H exp{W),4'-^'-^) 

= slp^^ZiM,;SZ). 

(4.16) 

Thus the one-7 expectation value is 

(7^)M.,^^ = 5^5%. (4.17) 



We next consider the two-7 function. We insert two 7s separated by t lattices 
in the direction along S . As same as above, the two-7 function is calculated 
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straightforwardly as 



aI a' 



= P-'' E ^M^MAfk' - J'l') + m - J'l))/P)) 

g 

(4.18) 
where W{T — t) and W{t) are the short expressions for the Ws for the intervals 
[t,T] and [0,t], respectively. We assumed for simphcity that the /i^-s at [0,t] and 
[t,T] are equal, and the H is defined as H{t;m;i,j) = Yl^=o^^Pi^^7k + ^27rA;j'/p). 

To evaluate (4.18) further, we will consider two cases separately. First consider 
the case k = k', I = I'. Then the Kronecker's deltas in (4.18) requires i = —i', j = 
—j', and we obtain 

p—l p—l 

= Z{Mg; 6Z)6l%^,6fl^,^,{ J2 exp{Th',,))~\ ^ exp{Thl,))-' 

r,s—0 r,s=0 

p—l P^l 

X ( J2 exp(tC, + (T - t)h'r,s-j))i Yl ^Mthi+i,s + {T- t)hl^s+j))- 

r,s—0 r,s—0 

(4.19) 
In the thermo-dynamical limit T ^ oo, we fix the separation t. Suppose h^^ ^^ 
and /i^^ „j have the largest real values among /i^^s and /i^ ^s (r, s = 0, ■ ■ ■ ,p — 1), 
respectively. Then we obtain in the thermo-dynamical limit 



■' J 



(4.20) 
Combining (4.20) and (4.17), one can see that the correlation lengths associated 
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to the 7^] iihj) 7^ (0)0)) ^^^ of order (5Z) '^ Then, in the hmit of TLFT, these 
correlation lengths diverge. 

Next consider the case {k,l) ^ {k',l'). Since then the Kronecker's deltas in 
(4.18) restrict i and i' to be zero, we obtain 

^(M,;7i^(t),7g;(0);5Z) = 5lp^^S^l5^lziM,;5Z). (4.21) 

Hence, 

(Tf^nU')''-'' = (7a^'^^(7g)^'^^. (4.22) 

This implies that the non-zero contributions to the connected two 7 functions 
of (A;, /) 7^ {k', I') will come only from the cluster deformations of the composite 
operators, and hence are of order (SZ) if the two 7s are separated by d lattices 
in the direction perpendicular to the S^. Thus the correlation lengths to that 
direction are of order —l/ln6Z, which vanish in the limit of the TLFT. 

We have obtained that, in the limit of the TLFT, the fluctuations associated to 
the operator 7 have the infinite correlation lengths in the direction along S^, but 
have zero correlation lengths in the directions perpendicular to S^. The fluctuations 
propagate only in the one-dimensional direction along the holes. We could not flnd 
any local fluctuations which propagate in all directions in Mg in the perturbed Zp 
model. 
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5. Summary, comments and discussions 

We investigated the thermo-dynamical natures near TLFTs, and found that 
they are in general on discrete phase transition points, and that they are on fixed 
points of renormahzation group transformations at least in a restricted sense. 

First we discussed the decomposition of a volume-dependent TLFT satisfy- 
ing a certain condition and showed that it is a direct sum of irreducible TLFTs 
with volume-dependent numerical factors. In the parameter space of the volume- 
dependent TLFTs we can discuss the phase structure and the renormahzation 
group fiow exactly, and showed that a TLFT is on a fixed point of the fiow and 
in general on a multiple first order phase transition point. The number of the 
different phases around a TLFT is equal to the number of the irreducible TLFTs 
of which the TLFT is a direct sum. 

To generalize the discussion out of the parameter space of the volume-dependent 
TLFTs, we introduced a kind of cluster expansion. By this perturbative scheme 
we found the same phase structure around a TLFT in the first order of the cluster 
expansion. Higher orders do not change the result qualitatively provided the base 
manifold is trivial. 

To specify the roles of the non-trivial topology of the base manifold and the 
physical states on non-trivial {D — l)-dimensional boundaries, we investigated the 
neighborhood of the Zp analogue of Turaev-Viro model in three dimensions. Then 
we found another phase structure, that is, the TLFT is on a higher order discrete 
phase transition point controlled by the physical states on boundaries with non- 
trivial topologies. 

We also studied the correlation lengths of fiuctuations associated to various 
boundaries in the cluster expansion near TLFTs. The correlation lengths of the 
fiuctuations associated to trivial local operators were shown to be zero in the limit 
of TLFT. On the other hand, in the Zp model, we found the fiuctuations whose 
correlation lengths diverge in the limit of TLFT. But these fiuctuations propagate 
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only along the topological defects. In this sense, a TLFT is on a point where only 
the fluctuations propagating along topological defects remain and have infinite 
correlation lengths. 

Comments are in order. We restricted our discussions in the cases that the 
numbers of the elements of the index sets Xi are finite. This restriction makes the 
numerical factors of the three-sphere functions Xi to be non-zero, and then we could 
obtain the decomposition theorem of a TLFT. But some continuum formulations of 
TFTs might allow A to be zero. The twisted N = 2 Landau-Ginzburg model '' with 
a certain super-potential is such a case. The TLFT corresponding to the model is 
discussed in the paper by Fukuma, Hosono and Kawai , where a certain infinite 
limit must be taken to show the correspondence. Moreover it was pointed out 
that the zero-coupling limit of 3D lattice QCD corresponds to the Ponzano-Regge 
model. Can we apply our arguments around the Ponzano-Regge model? It would 
be too stupid to argue that there are no local long-range fluctuations just around 
the zero-coupling limit of the 3D lattice QCD. Thus we insist that our discussions 
are applicable only to TLFTs with finite degrees of freedom, and new discussions 
are needed for the other types of TLFTs and TFTs. 

As we have shown, the physical states on S^^^ of a TLFT label the phases 
around a multiple first order phase transition point. The continuum correspon- 
dences of the physical states on S would be scalar observables of TFTs. Thus 
one could expect that such scalar observables would have special roles in the dy- 
namics of TFTs such as 'spontaneous symmetry breaking'. Such discussions were 
done by Kawamoto and Watabiki in their generalized Chern-Simons theory, 
which contains fields of forms of all orders. Here we cannot comment anything 
about it, since the applicability of our results is obscure as we mentioned in the 
last paragraph. Relating TLFTs to TFTs and developing analogous discussions 
for TFTs remain as future problems. 

The thermo-dynamical natures near a fixed point of a renormalization group 
transformation on a continuous phase transition surface can be well analyzed by 
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the renormalization group technique and renormahzable field theories. This is 
because of the universahty, or because some finite number of interacting long- 
range modes with some finite number of couplings control the system essentially. 
How about discrete phase transitions? In general we cannot expect universalities 
around such phase transition surfaces, so predictions are difficult. But when we 
analyzed the thermo-dynamical natures near TLFTs by the cluster expansion, 
the deviations could be projected to the physical states of TLFTs, and we could 
obtain some qualitative natures near TLFTs. In this sense, the thermo-dynamical 
natures near TLFTs are described by the physical observables of TLFTs. If there 
were a map to transform any point near a discrete phase transition surface to the 
neighborhood of a TLFT, one could talk about general discrete phase transitions. 
But the existence of such a map is obscure. Some difficulties of the renormalization 
group technique near a first order phase transition surface were pointed out . We 
can define a renormalization group fiow in the parameter space of the volume- 
dependent TLFTs, but will have to overcome difficulties in the outside. TLFTs 
can generate various systems with discrete phase transitions, but whether they are 
useful in the investigations of discrete phase transitions is not clear. 

Acknowledgements: We would like to thank H. Kawai, N. Kawamoto and A. Ya- 
mada for encouragement and the members of the elementary particle sections of 
Tokyo University for their hospitality. And we would like to thank H. Kawai and 
N. Kawamoto for reading this manuscript. The author is supported by the fel- 
lowship of the Japanese Society for the Promotion of Science for Japanese Junior 
Scientists, and this work is supported by the Grant-in-Aid for Scientific Research 
from the Ministry of Education No. 04-1324. 



45 



REFERENCES 

1. E. Witten, Commun. Math. Phys. 117 (1988), 353. 

2. M. Atiyah, Publ. Math. IHES 68 (1989), 175. 

3. E. Witten, Nud. Phys. B340 (1990), 281. 

4. E. Witten, Nucl. Phys. B311 (1988), 46. 

5. G. Moore and N. Seiberg, Phys. Lett. B220 (1989), 422. 

6. E. Witten, Nucl. Phys. B322 (1989), 629; E. Witten, Nucl. Phys. B330 
(1990), 285. 

7. E. Witten, Commun. Math. Phys. 121 (1989), 351. 

8. E. Witten, Commun. Math. Phys. 117 (1988), 353. 

9. J.F. Wheater, Phys. Lett. B223 (1989), 451. 

10. J.F. Wheater, Phys. Lett. B264 (1991), 161. 

11. T. Jonsson, Phys. Lett. B265 (1991), 141. 

12. T. Jonsson, Nucl. Phys.(Proc. Suppl.) B 25A (1992), 176. 

13. C. Bachas and P.M.S. Petropoulos, Commun. Math. Phys. 152 (1993), 191. 

14. M. Bordemann, T. Filk and C. Nowak, "Topological actions on 2-dimensional 
graphs and (graded) metrised algebras", preprint THEP 92/18. 

15. M. Fukuma, S. Hosono and H. Kawai, "Lattice topological field theory in 
two-dimensions", preprint CLNS 92/1173. 

16. B. Durhuus and T. Jonsson, "Classification and construction of unitary 
topological field theories in two dimensions" , preprint COPENHAGEN-MI- 
14-1993. 

17. A. Migdal, Sov. Phys. J.E.T.P. 42 (1975), 413. 

18. T. Eguchi and S.K. Yang, Mod. Phys. Lett A5 (1991), 1693. 

46 



19. G. Ponzano and T. Regge, in Spectroscopic and Group Theoretical Methods 
in Physics, ed. F. Bloch (North-Holland, Amsterdam, 1968). 

20. T. Regge, Nuovo Cimento 19 (1961), 558. 

21. V.G. Turaev and O.Y. Viro, Topology 31 (1992), 865. 

22. R. Dijkgraaf and E. Witten, Commun. Math. Phys. 129 (1990), 393. 

23. B. Durhuus, H. Jacobsen and R. Nest, Nucl. Phys.(Proc. Suppl.) 25A (1992), 
109; B. Durhuus, H. Jacobsen and R. Nest, Rev. Math. Phys. 5 (1993), 1. 

24. H. Ooguri and N. Sasakura, Mod. Phys. Lett A6 (1991), 3591. 

25. D.V. Boulatov, Mod. Phys. Lett A7 (1992), 1629. 

26. F. Archer and R.M. Wilhams, Phys. Lett. B273 (1991), 438. 

27. S. Mizoguchi and T. Tada, Phys. Rev. Lett. 68 (1992), 1795. 

28. F.J. Archer, Phys. Lett. B295 (1992), 199. 

29. G. Felder and O. Grandjean, "On combinatorial three-manifold invariants", 
ETH preprint (1992). 

30. S. Mizoguchi, Int. J. Mod. Phys. A8 (1993), 3909. 

31. H. Ooguri, Prog. Theor. Phys. 89 (1993), 1. 

32. V.G. Turaev, C.R. Acad. Sci. Paris 313 (1991), 395; V.G. Turaev, J. Diff. Geom. 36 
(1992), 35. 

33. H. Ooguri, Nucl. Phys. B382 (1992), 276. 

34. S. Chung, M. Fukuma and A. Shapere, "Structure of topological lattice field 
theories in three dimensions", preprint CLNS 93/1200. 

35. H. Ooguri, Mod. Phys. Lett A7 (1992), 2799. 

36. A.S. Schwarz, Lett. Math. Phys. 2 (1978), 247; G.T. Horowitz, Gom- 
mun. Math. Phys. 125 (1989), 417. 

37. N. Sasakura, Phys. Lett. B316 (1993), 329. 

47 



38. B.Ye. Rusakkov, Mod. Phys. Lett A5 (1990), 693. 

39. J.W. Alexander, Ann. Math. 31 (1930), 292. 

40. M. Gross and S. Varsted, Nucl. Phys. B378 (1992), 367. 

41. E. Verlinde, Nucl. Phys. B300 (1988), 360. 

42. D.V. Boulatov, V.A. Kazakov, I.K. Kostov and A.A. Migdal, Nucl. Phys. B275 
(1986), 641. 

43. See, for example, C. Itzykson and J.M. Drouffe, Statistical field theory, 
Cambridge (1989), vol.2, chapter 7. 

44. B. Nienhuis and M. Nauenberg, Phys. Rev. Lett. 35 (1975), 477. 

45. For review see, R. Dijkgraaf, H. Verlinde and E. Verlinde, "Notes on Topo- 
logical String Theory and 2D Quantum Gravity", Trieste Spring School 1990: 
91-156; 

R. Dijkgraaf, "Intersection Theory, Integrable Hierarchies and Topological 
Field Theory", NATO ASl: Cargese 1991: 95-158. 

46. N. Kawamoto and Y. Watabiki, Phys. Rev. D45 (1992), 605; N. Kawamoto 
and Y. Watabiki, Nucl. Phys. B396 (1993), 326. 

47. See, for example, A.C.D. van Enter, R. Fernandez and A.D. Sokal, Phys. Rev. Lett. 66 
(1991), 3253, and references therein. 



48 



FIGURE CAPTIONS 

1) Two lattices Mi and M2 with boundaries are glued at sub-boundaries b. 

2) A lattice M2 is obtained by identifying sub-boundaries b of Mi. 

3) Two lattices with dT as their simplically decomposed boundaries are glued 
at one of their T s of each boundary. This figure is for D = 3. 

The 2D lattice QCD is invariant under the flip move, and hence is an area- 
dependent TLFT. 

If the 6Zs come from the D-simplices shown in the figure, to estimate the 
contributions to the partition function, one has to use correlation functions 
of the TLFT on manifolds different from M. 

If the set Q is separated from the boundaries and the compliment Q near the 
boundaries, then W^ is zero. 

If the two boundaries are separated, the lower order terms of W^'^^ factorize. 

The (1,4) and (2,3) moves of the simplicial decomposition and their inverses. 
In three dimensions, all the simplicial decompositions of a manifold can be 
generated by these moves. 

9) The geometrical condition of t x Bg (t G S*^). 

10) r^^i is a physical state on the torus winding through the two neighboring 
holes k and I. 
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